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ABSTRACT 

Using the Submillimeter Array (SMA) we have imaged for the first time the 321.226 GHz, IO29 — 
ortho-H20 maser emission. This is also the first detection of this line in the Cepheus A high-mass 
star-forming region. The 22.235 GHz, 616 — 523 water masers were also observed with the Very Large 
Array 43 days following the SMA observations. Three of the nine detected submillimeter maser spots 
are associated with the centimeter masers spatially as well as kinematically, while there are 36 22 GHz 
maser spots without corresponding submillimeter masers. In the HW2 source, both the 321 GHz and 
22 GHz masers occur within the region of ~ 1" which includes the disk-jet system, but the position 
angles of the roughly linear structures traced by the masers indicate that the 321 GHz masers are 
along the jet while the 22 GHz masers are perpendicular to it. We interpret the submillimeter masers 
in Cepheus A to be tracing significantly hotter regions (600^2000 K) than the centimeter masers. 
Subject headings: ISM: individual (Cepheus A) — ISM: jets and outflows — masers — stars: formation 



1. INTRODUCTION 

Since its discovery in the interstellar medium nearly 
four decades ago (Cheung et al. 1969), the water maser 
emission at 22.235 GHz from the 616 — 523 transition 
has been extensively studied in a variety of sources, from 
circumstellar envelopes of late-type stars, star-forming 
regions and external galaxies (e.g., Elitzur 1992). This 
maser emission is understood to be arising from shocked 
gas with temperatures ~500 K and densities ~ 10 9 cm -3 
(Elitzur et al. 1989). 

Single-dish observations of the 22 GHz water maser 
emission have helped identify the early stages of massive 
star-formation (e.g., Sridharan et al. 2002; Brand et al. 
2003). Very Large Array (VLA) and Very Long Baseline 
Interferometry (VLBI) observations have allowed prob- 
ing of kinematics of gas in such regions with very high 
angular resolution of 0"08-0.5 mas (e.g., Claussen et al. 
1998; Furuya et al. 2000, 2005; Patel et al. 2000; Imai 
et al. 2000; Torrelles et al. 2001a, 2001b; Gallimore et 
al. 2003; Goddi et al. 2005; Goddi & Moscadelli 2006; 
Moscadelli et al. 2005). However, with observations of 
just the 22 GHz transition alone, it remains difficult to 
characterize the relative differences in the physical condi- 
tions in a given star-forming region. Nearly simultaneous 
observations of multiple masing transitions at compara- 
ble angular resolutions would be very helpful to allow 
probing of physical conditions of the masing gas such as 
its temperature and density, beyond just the lower limits 
required to satisfy the maser excitation. 

Except for the 22 GHz transition, and the 183.310 
GHz 3i3 — 220 transition, other transitions of H2O lie 
in the submillimeter band (Neufeld & Melnick 1991). Of 
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these, several submillimeter masers have been detected 
with single-dish observations from star-forming regions 
as well as late-type circumstellar envelopes (Menten et 
al. 1990a, 1990b; Melnick et al. 1993). Interferometric 
observations of some of these transitions are only now 
possible with the availability of the Submillimeter Array 
(SMA; Ho et al. 2004, Humphreys et al. 2005). 5 

Cepheus A is a well-studied high-mass star-forming re- 
gion in the Cepheus OB3 complex of molecular clouds 
(Sargent 1977). Several compact objects in this region 
have been identified from radio continuum observations 
at centimeter wavelengths with the VLA (hereafter de- 
scribed as HW sources: Hughes & Wouterlout 1984; 
Garay et al. 1996; Martin-Pintado et al. 2005). The 
22 GHz water masers have been observed towards sev- 
eral of these HW sources (Torrelles et al. 1996, 2001a, 
2001b; Vlemmings et al. 2006) but observations of sub- 
millimeter masers have the potential of allowing a study 
of differences in the physical conditions in these sources, 
and hence may lead to a better understanding of their 
evolutionary differences. 

Here we report the first detection of the submillime- 
ter wavelength water maser emission in the Cepheus-A 
region. The fact that this emission is detected at all, sug- 
gests the presence of highly excited gas due to the very 
high energy levels of the submillimeter maser transitions. 

2. OBSERVATIONS 

The Cepheus A star-forming region was observed with 
the SMA on 2004 August 30 using seven antennas in 
the extended configuration with a maximum baseline 
length of 220 m. The phase center was a(2000) = 
22 ft 56 m 17.971 s , 5(2000) = +62°01'49."279. We used a 
tuning of 321.226 GHz to center the water maser lines 
in the lower sideband (LSB). The quasars BL Lac and 
0102+584 were observed for 5 minutes each between ev- 
ery cycle of 20 minutes on the main source. The spectral 
band-pass was calibrated using observations of Uranus 

5 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics, and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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Fig. 1. — VLA map of 1.3 cm continuum emission from HW2 
and HW3c (inset). The positions of sub-millimeter water mascrs 
observed with the SMA are shown with error bars representing 
the formal uncertainties. The position offsets were obtained by 
fitting two-dimensional elliptical Gaussians to integrated intensity 
maps over channels in which the maser emission appeared (using 
the IMFIT task in Miriad). The position offsets are with respect 
to the absolute coordinates: a(2000) = 22 h 56 m 17.971 s , <5(2000) = 
+62°01'49."279. The crosses show the positions of the 22 GHz 
water masers observed with the VLA (this Letter). The errors in 
the positions of the 22 GHz water masers are smaller than the size 
of the symbols shown. Only three of these centimeter wavelength 
masers appear to be correlated with the submillimeter masers in 
both position and velocity (see Table 1). Only HW2 and HW3c 
show the presence of submillimeter water masers among all the 
compact continuum sources in the entire Ccphcus-A star-forming 
region (over ±5"). The VLA beam is 0."09 X 0."07 (P.A. 14°). 
The SMA beam is 0."8 X 0."7 with a P.A. of -80°. 



and Saturn and absolute flux calibration was done us- 
ing Neptune and the quasar 2232+117. Results from 
the upper sideband data (centered at 331 GHz, includ- 
ing CH 3 CN lines and continuum emission) have been re- 
ported earlier by Patel et al. (2005). Weather conditions 
were excellent during the observations with relative hu- 
midity ~ 10% and t 2 25GHz « 0.08 {t 32 ighz ~ 0.69), 
measured at the nearby Caltech Submillimeter Observa- 
tory. The track was ~10 hr long with on-source integra- 
tion time of ~6.5 hr. T sySt DSB varied from 320 to 500 
K. The visibility data were calibrated using the Owens 
Valley Radio Observatory's MIR package and imaging 
was done using the Miriad package. We applied the 
self-calibration solutions obtained from the continuum 
emission (of ~ 2 Jy beam -1 ) from the HW2 source ob- 
tained from the LSB data, to the spectral-line data of the 
submillimeter wavelength water masers. With uniform 
weighting, the synthesized beam was 0."8 x 0."7 with a 
position angle (P. A.) of -80°. We expect the maximum 
error in our absolute astrometry to be ~ 0."1, based on 
positions of quasars mapped in our SMA observations. 
We estimate an uncertainty of ^20% in the absolute flux 
scale in the SMA data. 
We also observed the 22.235 GHz water masers in 
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Fig. 2. — Spectra of submillimeter and centimeter wavelength 
water maser emission from Cepheus A HW2 and HW3c. 

Cepheus-A with the VLA of the NRAO 6 in the A con- 
figuration on 2004 October 12. These observations were 
made in snapshot mode with a total integration time of 
about 15 minutes on-source. The rms noise level in our 
VLA observations is about 3.3 x 10~ 3 Jy beam -1 in chan- 
nels free of maser emission and ~4.5xl0 -1 Jy beam -1 
in the channel with the strongest maser spot. The line 
and continuum data at 1.3 cm were obtained simultane- 
ously and the astrometric registration between the maser 
and continuum emission is better than 0."01. These data 
were reduced using the NRAO AIPS package. 

3. RESULTS 

Maps of the 1.3 cm continuum emission observed with 
the VLA (on 2004 October 12), with an overlay of the 
22 GHz water masers and the 321 GHz water masers 
observed with the SMA are shown in Figure 1. The 
321 GHz water masers were detected only in HW2 and 
HW3c. The HW2 thermal jet has been observed earlier 
in several epochs from 1991 to 2004 and proper motions 
of ~ 500 km s _1 of this jet were recently s reported by 
Curiel et al. (2006). The position angle of the jet of ~ 
45° is nearly perpendicular to the proposed circumstellar 
disk seen as a flattened elongated structure in continuum 
emission at 900 ^m and CH 3 CN line emission (Patel et al. 
2005). A large majority of the 22 GHz masers are roughly 
along the direction of the disk (which is a larger struc- 
ture of 1" beyond the scale of Fig. 1), consistent with 

6 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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that found earlier by Torrelles et al. (1996, 2001a,2001b) 
in their VLA and Very Long Baseline Array (VLBA) ob- 
servations. The 321 GHz masers in HW2 appear to lie 
along the outflow jet. A total of 39 masers were detected 
at 22 GHz, only those associated with HW2 and HW3c 
are shown in Fig. 1. 

Figure 2 shows the spectra of 22 and 321 GHz masers in 
HW2 and HW3c. With respect to the systemic velocity 
of the ambient molecular cloud in Cepheus A (Vlsr — 
-11 km s _1 ; Gomez et al. 1999), the strongest emission 
in submillimcter masers associated with HW2 appears 
to occur at relatively higher velocities, whereas the 22 
GHz masers are significantly weaker at these velocities. 
The brightest 22 GHz masers have flux density of ~1400 
Jy, while the brightest 321 GHz maser is ^4 Jy. The 
22 GHz masers are about a factor of 1000 times brighter 
than the 321 GHz masers in general. The observed char- 
acteristics of the submillimcter masers are listed in Table 
1 along with corresponding values for the three spatially 
and kinematically associated 22 GHz masers. 

4. DISCUSSION 

From the sizes of the 321 GHz emission spots in channel 
maps and the observed flux densities, we estimate a lower 
limit in the brightness temperature to be ~200 K. Since 
we do not have a reliable estimate of the sizes of these 
spots, we are unable to check whether the emission is 
due to maser action based on the brightness temperature. 
The narrow line-widths of the 321 GHz emission (0.3-1 
km s _1 ) compared to the widths of the CH 3 CN lines, 
strongly suggest that the 321 GHz lines arc indeed due 
to maser emission. 

The first astronomical detection of the IO29 — 936 tran- 
sition of H2O masers was made by Menten et al. (1990a) 
with the Caltech Submillimeter Observatory 10.4 m tele- 
scope. Shortly after the discovery of the 321 GHz masers, 
Neufcld & Melnick (1990, hereafter, NM90), showed that 
the observed 22 and 321 GHz masers can be explained 
by collisional excitation within the same volume of gas. 
NM90 presented models of the emissivity ratio R of 22 
GHz / 321 GHz maser luminosities as a function of £ and 
T, where £ is a combination of parameters involving the 
gas density, abundance of water and the magnitude of 
velocity gradient in the slab of gas that is masing (see 
equation 1 of NM90). In a subsequent paper, Neufeld & 
Melnick (1991) presented more detailed theoretical calcu- 
lations including 349 rotational states of ortho and para 
water under a range of £ and radiation field. The 321 
GHz transition is among the 20 other theoretically pre- 
dicted maser transitions in the frequency range of 183- 
1542 GHz. Many of these transitions are shown to be 
approaching in brightness to the 22 GHz maser, under 
certain ranges of T and £. Yates et al. (1997) carried out 
similar calculations for several submillimcter H2O maser 
transitions, including the 321 GHz line. Their range of 
physical conditions is wider than in NM90 and Neufeld 
& Melnick (1991), and they also considered the effects of 
radiation from dust grains. According to these models, 
the 321 GHz maser is strongly inverted under more re- 
stricted conditions than the 22 GHz line (see Fig. 3 of 
NM90 and Fig. 4a of Yates et al. 1997). In fact, strong 
inversion of the 321 GHz transition requires T^ > 1000K. 
Assuming a temperature of 1400 K, the greatest gain in 
the 321 GHz line occurs for gas density n(H 2 ) in the 



range of 4-6xl0 8 cm~ 3 (Yates et al. 1997). 

All of these model calculations have a central assump- 
tion that the submillimeter and centimeter wavelength 
masing transitions arise in the same volume of gas. This 
implies that the 22 GHz masers have a broader range of 
physical conditions in temperature and densities (Yates 
et al. 1997). (There can be 22 GHz emission and no 
accompanying 321 GHz masers in the same region but 
not vice- versa, according to these models). We can con- 
front these models with our maps of 22 and 321 GHz 
masers, to first ask if these masers arise in the same vol- 
ume of gas, and if they do, what are the implications of 
the observed luminosity ratio. 

Figure 1 and Table 1 show that we have two out of 
eight masers in HW2 that appear to be coincident in 
both spatial and velocity neighborhood (< 0."1 and < 3 
km s _1 ) with 22 GHz masers. In addition, we have one 
more region in the HW3c which is located ^3" south 
of HW2, that appears to have masers at both 22 GHz 
and 321 GHz arising from the same region. If we use a 
more stringent criteria for coincidence of masers in line- 
of-sight velocity to be within 1 km s , based on the 
typical line widths of the maser emission, we would have 
only one maser spot (occurring in HW3c; number 9 in 
Table 1) that would be coincident. The coincidence of 
masers 1 and 7 in HW2, might therefore be considered to 
be tentative, given the larger difference in their velocities. 
We speculate that this difference in velocity is related to 
the fact that HW2 is a dynamically more active region 
compared to HW3c. 

The emissivity ratio R (22/321 GHz) for these masers 
listed as numbers 1, 7, and 9 in Table 1 are 9.2, 0.02, 
and 0.18. The value of R is roughly expected to decrease 
with temperature, given that E/k — 643K for the 22 
GHz transition and 1861 K for the 321 GHz line. In this 
way, a much hotter gas is required for the submillimc- 
ter wavelength maser excitation. Assuming that both 22 
and 321 GHz masers are saturated, the luminosity ratio 
is then the same as the ratio of rate coefficients as defined 
by NM90 and plotted in their Fig. 3. We find that for 
maser- 1 (Table 1), the observed ratio implies the kinetic 
temperature to be in the range of 500-2000K, depending 
on the value of £. According to NM90, the typical value 
of log £ would be in the range of -1.25 to -0.25 for the 
shocked gas in star-forming regions and their Figure 3 
would then imply a temperature of ~ 600K for the re- 
gion where maser 1 occurs. Masers 7 and 9 have much 
smaller values of R, much less than 1. Even for R = 5, 
according to NM90, such masers would presumably ex- 
ist behind slower non-dissociative shocks in which the 
kinetic temperature of the gas is ~ 1000K. Again, ac- 
cording to figure 3 of NM90, this value of R < 1 implies 
gas temperatures as high as 2000K for a wide range of £. 

With respect to the remaining six masers (listed as 
numbers 2-6 and 8 in Table 1) that show emission only 
in the 321 GHz and not in the 22 GHz transition, we 
cannot apply any of the existing theoretical models to 
these R — > cases that imply that the 22 GHz masers 
are quenched but the 321 GHz emission remains strong. 
According to Figure 3 of NM90, which has the lowest con- 
tour corresponding to R = 1.5, the regions with R — > 
may imply temperatures greater than 2000K. On the 
other hand, we cannot completely rule out the effect of 
time variability for the absence of 22 GHz masers corre- 
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sponding to some the observed 321 GHz masers. 

The 22 GHz masers are well known to be time- variable. 
Previous VLBA observations of Cepheus A (Torrelles et 
al. 2001a, 2001b) and IRAS 21391+5802 (Patel et al. 
2000), show that several masing spots disappear over 
timescales of 1 month. These masers however, tend to 
be relatively weaker. The theoretical models as in NM90 
predict strong emission at 22 GHz corresponding to the 
321 GHz masers. If these corresponding 22 GHz masers 
existed during the epoch of the submillimeter observa- 
tions, they are unlikely to have been completely extin- 
guished about a month later in the centimeter wave- 
length observations. 

Much less is known about the variability of the 321 
GHz masers. In late-type stars, near simultaneous ob- 
servations of 22, 321, and 325 GHz masers were carried 
out by Yates & Cohen (1996) over four epochs separated 
by ~ 20 days. They find the 321 GHz masers to be the 
most variable, in terms of both changes in brightness - 
by a factor of 2 - 10 and the shortness of time scale, by 
~20 days. However, it remains unclear how far we can 
carry forward these findings to the case of water masers 
associated with star-forming regions. 

The 22 GHz masers that are not associated with the 
321 GHz masers are likely to be arising in relatively 
cooler regions among the various HW sources in Cepheus 
A, compared to HW2 and HW3c. The fact that the sub- 
millimeter wavelength masers in HW2 are found along 
the major axis of the jet, lead us to propose that they 
are arising in hotter gas regions due to the impact of the 
jet. 

5. CONCLUSIONS 

We report the first detection of 321 GHz water masers 
in the high-mass star-forming region Cepheus A. The 
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SMA observations at ~ 0."8 have sufficient spatial reso- 
lution to allow a study of the association of these masers 
with the multiple radio continuum sources such as HW2 
and HW3c (separated by about 3"), and also with the 
disk-outflow system of HW2 which has a size scale of 
about 1" (725 AU). We had near-simultaneous 22 GHz 
water maser observations made with the VLA that show, 
in general, a poor agreement between the locations of 
submillimeter and centimeter water masers: only three 
out of nine masers agree in both position and line-of-sight 
velocities. Moreover, the submillimeter masers are found 
to lie along the outflow jet whereas the 22 GHz masers 
are primarily found perpendicular to it. The mapping of 
submillimeter wavelength masers with the SMA, along 
with near-simultaneous mapping of the 22 GHz water 
masers with the VLA, opens the possibility of studying 
the physical conditions of density and temperature of the 
gas associated with high-mass star-forming sites (where 
conditions are likely to be suitable for excitation of sub- 
millimeter water masers), at very high spatial resolution 
(~ 100 AU). 
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TABLE 1 

Summary of detected masers 



321 GHz MASERS a 22 GHz MASERS b 



No. 


Aa 


A<5 


Flux 




Aa 


A<5 


Flux 










Density 


Vlsr 






Density 


Vlsr 




// 


// 


(Jy) 


(km s -1 ) 


// 


// 


(Jy) 


(km s _1 ) 


1 


0.096 


0.359 


2.8 


-19.0 


0.068 


0.378 


25.9 


-21.1 


2 


0.012 


0.338 


4.8 


-21.0 










3 


-0.037 


0.276 


2.9 


-11.0 










4 


-0.016 


0.193 


2.7 


-4.5 










5 


0.005 


0.182 


3.0 


-1.0 










6 


-0.023 


0.180 


6.5 


-3.5 










7 


-0.030 


0.168 


5.1 


1.2 


-0.004 


0.167 


0.09 


4.0 


8 


-0.023 


0.136 


1.7 


-3.8 










9 C 


-0.115 


-3.131 


3.7 


-5.5 


-0.005 


-3.133 


0.67 


-5.6 



a Typical uncertainty in SMA flux values is 0.4 Jy. Formal uncertainty in Gaussian fitted positions is ~ 30 mas (see error-bars on position 
offsets showned in Fig. l). b The 22 GHz masers listed here are a subset of all the masers detected in the VLA observations. Only those 
masers which are most likely to be associated with the 321 GHz masers (closest occurrences in position and velocity) are tabulated here. 
Typical uncertainty in VLA flux values is 10 mjy and formal uncertainty in positions is ~ 1 mas. c This maser is associated with HW3c, 
all others are associated with HW2. 



